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Early Dark Energy Cosmologies 
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We propose a novel parameterization of the dark energy density. It is particularly well suited 
to describe a non-negligible contribution of dark energy at early times and contains only three 
parameters, which are all physically meaningful: the fractional dark energy density today, the 
equation of state today and the fractional dark energy density at early times. As we parameterize 
fid (a) directly instead of the equation of state, we can give analytic expressions for the Hubble 
parameter, the conformal horizon today and at last scattering, the sound horizon at last scattering, 
the acoustic scale as well as the luminosity distance. For an equation of state today wq < — 1, 
our model crosses the cosmological constant boundary. We perform numerical studies to constrain 
the parameters of our model by using Cosmic Microwave Background, Large Scale Structure and 
Supernovae la data. At 95% confidence, we find that the fractional dark energy density at early 
times fi d < 0.06. This bound tightens considerably to fi d < 0.04 when the latest Boomerang data 
is included. We find that both the gold sample of Riess et. al. and the SNLS data of Astier et. al. 
when combined with CMB and LSS data mildly prefer wo < — 1, but are well compatible with a 
cosmological constant. 

PACS numbers: 98.80.-k 



INTRODUCTION 

Current observations 0, 0H II E| favor some form 
of dark energy that today comprises roughly 70% of the 
energy density of our Universe. One fundamental issue 
is whether dark energy is a true cosmological constant 
or time evolving In recent years, the notion of 

an evolving dark ene rgy has been cast in various param- 
eterizations El El El El El El E3 of the equation 
of state w(a) = p/p of dark energy. Yet such parame- 
terizations are ill-suited to catch an intriguing possible 
feature of dark energy, namely that it could be present 
at an observable level even from as early as Big Bang 
Nucleosynthesis on. Such models would leave their im- 
prints on the Cosmic Microwave Background , cosmic 
structure Il8l 1 191 and maybe even Big Bang Nucleosyn- 
thesis |2CJ, |2l[ |22j • There are some good points in favor 
of this scenario: if one links the presently small energy 
density of dark energy to the age of the Universe, one 
is led to attractor solutions E|- Such scenarios occur 
in attempts to solve the cosmological constant problem 
from the point of view of dilatation symmetry and 
also in string theories. 

Instead of parameterizing w(a), we parameterize fid (a) 
directly. This will prove advantageous for two reasons: 
firstly, the amount of dark energy at early times is then a 
natural parameter and not inferred by integrating w(a) 
over the entire evolution. Secondly, since the Hubble 
parameter is given by 

g 2 (a) = ^ na - 3 + ^ cl q- 4 
Hi l-fid(a) ' (> 

a simple, analytic expression for f2d(a) enables us to com- 
pute many astrophysical quantities analytically. In the 
above, f2° cl is the fractional energy density of relativis- 
tic neutrinos and photons today, f2^ is the matter (dark 




FIG. 1: Evolution of the fractional dark energy density fid(-z) 
and the equation of state w(z) as a function of redshift. The 
solid (black) curve depicts the behavior for a ACDM cosmo- 
logical constant model, in which u>o = — 1 by definition and 
the amount of dark energy at early times tends to zero. In 
contrast, the dashed (blue) and dotted (red) curves are early 
dark energy models described by our parameterization ©. 
For the models depicted, we chose wo = — 1 and fi d = 0.01 
(blue, dashed), fi d = 0.07 (red, dotted) respectively. In ad- 
dition, we plot the equation of state w for fi d = 0.01 (blue, 
dashed-dotted) and fijj = 0.07 (red, dashed-double-dotted) 



and baryonic) fractional energy density and we assumed 
a flat Universe. For any parameterization of f2d(a), the 
equation of state can of course be inferred from f2d(a) via 
an analytic relation (see Equation JSJ). 
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PARAMETERIZING DARK ENERGY 

As said, we would like to parameterize fid (a) to catch 
the important feature of early dark energy. In addition, 
the parameterization must involve only a restricted num- 
ber of - physically meaningful - parameters (we will only 
require three: fi^, wq and fi^, see below). Equation © 
is our proposal which is derived from the following sim- 
ple considerations: we start by observing that neglecting 
radiation, the fractional energy density of a cosmological 
constant evolves as 
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fiH, 
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Here, a is the scale factor normalized to a = 1 today, fij 
and fi^ are the fractional densities of dark energy and 
dark matter today and we assume a flat universe, i.e. 
fi^ + fijj = 1. The first generalization of this formula is 
to allow w ^ — 1, which is achieved by 0] 
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FIG. 2: Constraints on the parameters wo and fi§ for different 
combinations of data sets. The blue region corresponds to 
the SNe la compilation of Riess et. al., the green region to 
WMAP + VSA + CBI + SDSS. The constraints obtained 
when combining all of these sets are shown in yellow. The 
result of adding the Boomerang data to this combined set 
is depicted in red. Black (white) lines enclose 95% (68%) 
confidence regions. 



A straightforward attempt to include early dark energy is 
to simply add a term that is switched on at high redshifts 
and gives a basic contribution of the requested level 



fid (a) 



fi° 



fi5 + fi° a 3w ° 



«d(l-a Q )- 
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where a > is a parameter. It turns out, however, that 
this is insufficient, because the evolution of fid is con- 
nected to the equation of state w by the relation |2^| 



3w 



fi d (l - fid) = -dfid/dlna, (5) 



where a eq is the scale factor at matter-radiation equality. 
Demanding that w(a = 1) = wq, i.e. that the parameter 
wq should indeed have its usual meaning, one is led to 
conclude that an additional term is needed in the numer- 
ator and that a — —3wq- This yields the final form of 
our parameterization, namely 



n i \ d "d 
fid(a) = 
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In terms of the equation of state, going from today to 
the past, w starts at wq. It then crosses over toio«0 
during the matter dominated era. Defining the cross- 
over as w = wq/2 and using Equations JSJ and ijfjj and 
working to leading order in fi^/fi[j, one obtains the cross- 
over scale factor 



fi° 
Si d 



[1 - fi°]fi§ 



(7) 



We see from Equation {7|), that increasing fijj increases 
the cross-over scale factor, as expected. Likewise, a more 
negative wq also increases a c , i.e. cross-over of w occurs 
more recently. Finally, in the radiation dominated era, 
w tends to 1/3. 

It is worth mentioning that our parameterization © 
is a monotonic function of a for as long as 



fi" 

O e < d 



(8) 



i.e. it stays monotonic even for rather large fi^. Hence, 
as w = during matter domination, we see that 
for wq < —l, w(a) will cross the cosmological constant 
boundary [13, HH EE E3, E3| . The evolution of fi d and 
the corresponding w are depicted in Figure ^ 



ANALYTIC RESULTS 

A direct parameterization of fid (a) removes one in- 
tegration otherwise necessary to compute the luminosity 
distance, sound horizon etc., because the Hubble parame- 
ter H is given by Equation Q . Many cosmological quan- 
tities can then be computed analytically; the luminosity 
distance d^, the horizon today and at last scattering To 
and T/ s , the sound horizon r s and the acoustic scale Ia- 

To derive an expression for the horizon today To, we 
neglect radiation in Equation Q , which leads to an error 
of less than 1%. Combining Friedmann's Equation today 
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and at arbitrary scale factor yields 



da 



— ) =H$(n° m a + n d (a)) 



(9) 



Inverting, drawing the root and separating the variables, 
one gets 



dr = Hn 1 / da 



l-fi d (a) 



(10) 



Using our parameterization Equation JSJ and substitut- 
ing y = a -3 ™ , one obtains 
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Expanding the root in the numerator yields 

_i fif 1 ^ fi-f [i-y] 2 

TO 



3wo-ffo\/^m JO 



(11) 



dy. 



We then split the integral in leading and next to leading 
order 
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to 
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Both integrals yield hypergeometric functions 2-F1 (up to 
r functions) . The final result may hence be expressed as 
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In the limit that fl^ — ► 0, we recover the result for con- 
stant w (lif . The luminosity distance is computed in 
quite the same manner. Neglecting radiation, it is given 
by 

7 . . . . f* dz' . . f 1 a- 2 da 
d L (z) = (1 + z) =(l + z ) 

-1 



/o H(z) 
a-V 2 da 



1 - fi§[l - a- 3u, o] 2 



(15) 



The integral is in fact identical to the one in Ullfl. but 
with a(z) instead of as lower boundary. It seems feasi- 
ble to derive an expression for cIl(z) in terms of hyper- 
geometric functions as well 30]. Yet such an expression 



would be rather lengthy. As the evaluation necessitates 
numerical methods in any case, we leave Equation 1|15|) 
as it is. 

At high redshift, fi^ simply scales the Hubble parame- 
ter by a constant factor. The Friedmann Equation (this 
time inclu ding radiation) then yields the horizon at last 
scattering 
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Here, ai s is the scale factor at recombination. 
The sound horizon is given by 

dr 



r s (a) = / da—c s , 
Jo da 



(17) 



where the speed of sound is c s 



3(1 + R~ 



Here 

R = |— is the photon to baryon ratio. At redshifts 
where the speed of sound is appreciable and r s (a) receives 
contributions, Qd(a) ~ fi d holds very well and we get 
again from the Friedmann Equation that 



da 



Ho J ^/Zin^a + n^il + R- 1 )' 
which can be integrated (similarly to |3lJ) to 
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Here, Ri B and R equ . is the photon to baryon ratio as de- 
fined above at last scattering and matter-radiation equal- 
ity respectively. 

Finally, with r s , tq and Tu given, the acoustic scale 



TO - Tls 



(20) 



lA = 7T 

can easily be computed. 



CONNECTION TO EXISTING w{a) 
PARAMETERIZATIONS 



It may be worthwhile to relate our parameterization to 
the w(a) parameterization of (l2^ 



w Coras -(a) =w + (w m - wo) x 



1 + exp(a c /A) 
1 + cxp([a c — a]/ A) 
1 -exp([l -a}/ A) 



1 -exp(l/A) 



(21) 
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FIG. 3: The logarithm of the potential, ln(U) as a function 
of ip in units of Mp (solid line). The exponential potentials 
during radiation and matter domination are indicated by the 
dashed (blue) and dotted (red) line. The exponent in these 
cases is given by Equation 1231 . For this plot, we used — 
0.03 and matter-radiation equality is at ip equ . = —2.66. In the 
recent universe, the potential flattens, leading to w — > —0.99 
which we picked for this plot. 

This versatile parameterization is characterized by the 
equation of state of dark energy today Wq, the dark en- 
ergy equation of state during earlier epochs w m , a cross- 
over scale factor a c and a parameter A controlling the 
rapidity of this cross-over. In our case, wq has the same 
meaning, and w m — always for our parameterization. 
We might roughly identify a c in (|2*T)l with Q). Using 
equations © and © and working to leading order in 
ri|/il°, we also get an estimate of the width of the tran- 
sition 

A c r» (l - 3^) a c . (22) 

We caution the reader that the accuracy of these expres- 
sions varies and hence the relation of our parameteriza- 
tion to that given in Equation (|21|l should be seen as a 
semi-quantitative statement. 

THE CASE OF SCALAR DARK ENERGY 

If one uses our parameterization to describe the evo- 
lution of an (effective scalar dark energy field, then 
the constant Sid ~ ^d at high redshifts implies an expo- 
nential potential for the scalar field. Indeed, the value of 
Sl§ is that of an attractor solution for the potential Q 

V{(p) = M$ exp (-^3[l + w backg ]/n e d i^j , (23) 

where Mp is the reduced Planck mass and Wbackg is the 
equation of state of the components other than dark en- 
ergy, i.e. Wbackg = 1/3 during radiation domination and 
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FIG. 4: Upper left panel: Likelihood distribution for the equa- 
tion of state wo today for a phantom model with = 1. The 
solid line is for WMAP + VSA + CBI + BOOMERANG + 
SDSS + SNe la, the dashed (blue) line is without Sne la data. 
At lcr confidence, the full data gives wo = —1.08 + 0.14 — 0.17. 
Upper right panel: Likelihood distribution for the scalar spec- 
tral index n s for a scalar field model. The solid line is for 
WMAP + VSA + CBI + BOOMERANG + SDSS + SNe 
la, the dotted (red) line is without Boomerang. At la confi- 
dence, n B = 0.98 + 0.04 — 0.03. Lower left panel: Likelihood 
distribution for for a scalar field model. The solid line is 
for WMAP + VSA + CBI + BOOMERANG + SDSS+SNe 
la, the dotted (red) line is without Boomerang. At 2a con- 
fidence, < 0.04 and Sl^ < 0.06 respectively. Lower right 
panel: Comparison of the likelihood distribution for the equa- 
tion of state wo today for the gold set from Riess et. al. only 
(green, dot-dashed line) and the SNLS data only (blue, dotted 
line), and their combination with the CMB+LSS data: the 
dashed line (red) corresponds to CMB+LSS+Riess et. al., 
the solid (black) line to CMB+LSS+SNLS. 

Wbackg = during matter domination. As Sl| = const. 
in our parameterization, we conclude that in terms of a 
scalar field potential, we have four "phases". The po- 
tential (1) is an exponential potential during radiation 
domination, (2) moves over to a slightly less steep ex- 
ponential potential during recombination, (3) stays with 
that exponential potential during matter domination and 
at late times (4) the potential flattens or - equivalently 
- the kinetic term |33J, |3J] changes. In Figure |3 we ex- 
emplify this behavior by plotting ip in units of Mp vs. 
ln(U). For this plot, we have normalized ip such that 
today <p — 0. 

As said, our parameterization extends to phantom 
crossing models. For wq < —1, there is a crossing from 
the phantom to the "canonical" dark energy regime. In 
this case, however, it seems impossible to describe the 
behavior using a single scalar field [H |H El IH 113 • 
For simulations of phantom crossing models, we chose to 
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fix the rest frame speed of sound of dark energy to the 
speed of light c 2 = 1, which is compatible with the "mi- 
croscopic" behavior of scalar dark energy models. Yet, 
from the point of view of fundamental physics, it seems 
ill-suited, because the entropy generation rate V diverges. 
We nevertheless compute CMB and LSS constraints for 
phantom crossing models using c 2 — 1, because the de- 
tailed choice of c 2 seems to be of no significance for phe- 
nomenological studies using current data 28]. 

SIMULATION RESULTS 

We added the parameterization © to Cmbeasy |8f| 
and computed constraints on cosmological parameters 
using a Markov Chain Monte Carlo approach [3(| with 
the AnalyzeThis! package [13 • The parameters we 
chose are the matter energy fraction ^l^h 2 , the baryon 
energy fraction f^ft 2 , the hubble parameter h, optical 
depth r, scalar spectral index n Sl initial scalar ampli- 
tude A s using the observationally relevant combination 
ln(10 10 ^4 s ) — 2r as well as the dark energy parameters 
uuq and fl^. We chose flat priors on all parameters and 
ran two kinds of models separately: scalar field dark en- 
ergy models with wq > — 1 and phantom crossing models 
with wo € [—5,0] and speed of sound c 2 — 1. We com- 
pare the predictions of the models with several combi- 
nations of data sets. The first set, CMB+LSS, consists 
of CMB data from WMAP 0L VSA [| and CBI Q, as 
well as LSS data from SDSS [fL For SNe la, we use the 
compilation by Riess et. al., |2(. Motivated by recent 
comparisons of different SN la data sets 0, [3!| , we also 
take the first-year data from the SNLS [l| into account 
when considering phantom crossing models. 

Turning first to the case of scalar field dark energy, 
our main results are shown in Figure Eland Table [i] We 
see that the SNe la data and the CMB+LSS data give 
orthogonal information. While supernovae are insensitive 
to the amount of early dark energy but do constrain 
wo, the opposite is true for the CMB+LSS set. Figure[21 
also shows the likelihood contours obtained when adding 
the data from the 2003 flight of Boomerang [iol lil lij] . 

The stronger constraints on the amount of early dark 
energy from Boomerang are due to the particular sup- 
pression of power caused by fi^: the presence of dark 
energy suppresses the growth of linear fluctuations that 
are inside the horizon. The smaller the scale, the longer 
it suffered the suppression. This leads to a scale depen- 
dent red tilt of the spectra for all modes that enter after 
the scale of equality k < k equ , . All modes k > k equ , were 
inside the horizon before equality and are suppressed by 
the same factor 0]. An increase of is hence par- 
tially degenerate with a decrease of n s . Now Boome rang 
alone prefers a rather red spectral index n s = 0.86 [44j, 
whereas WMAP tends more towards n s ~ 1 (see also the 
upper right panel of Figure . As far as early dark en- 



ergy models with appreciable fijj are concerned, WMAP 
forces them to have n s close to 1, at least slightly larger 
than a comparable model with vanishing fl^. The rela- 
tive lack of power of Boomerang at high multipoles does 
then disfavor those models with large f2J and rather blue 
n s . 

Considering phantom crossing models, we use the 
CMB+LSS data, and the SNe la data from either Riess 
et. al. or the SNLS. In Table ITT1 we summarize the con- 
straints. The comparison of the results for the equation 
of state parameter today, wq, are shown in the upper left 
and lower right panels of Figure^] We find that the SNLS 
data alone, in contrast to the set from Riess et. al., do 
not favor a value of wq < 1, which agrees with the find- 
ings in [23. This situation changes when combining the 
supernovae data with CMB+LSS, since the CMB essen- 
tially fixes ^m^ 2 - Without a free both supernovae 
sets prefer a w of slightly less than -1 [Ifl Eft I47L Esj . 

CONCLUSIONS 

We introduced a new parameterization of dark energy. 
It is particularly well suited to describe models in which 
the dark energy density is non-negligible in the early Uni- 
verse. Working with £!d(a) instead of w(a) facilitates the 
computation of quantities such as horizons, the luminos- 
ity distance and the acoustic scale. Using Markov Chain 
Monte Carlo simulations, we constrained the parameters 
of our model using the latest CMB, Sne la and LSS data. 
While the CMB is more sensitive to the opposite is 
true for Sne la, which are more sensitive to wq. Adding 
the recent high multipole Boomerang data tightened the 
upper bound on fi^ considerably to £1^ < 0.04. However, 
our analysis only included linear fluctuations. As early 
dark energy models lead to more non-linear structure at 
higher redshifts compared to a cosmological constant > 
it might well be that the upper bound turns into a de- 
tection in the future. 

Acknowledgments We would like to thank Ben Gold 
for his helpful comments on the preprint. 

APPENDIX 

The upper bound for f2J with our parameterization in 
the scalar dark energy case is considerably higher than 
the one found in 0] . This is a consequence of the longer 
period of time for which our parameterization mimics a 
cosmological constant, so that the influence of fijj on the 
late-time evolution of the Universe induced by the choice 
of parameterization is minimized. It is then natural to 
ask if the bounds on fJ^j in our model might relax if f2d(a) 
behaved like a cosmological constant for higher redshifts. 
To answer this question, we generalize our parameteriza- 
tion Equation JfJJl by introducing an additional parameter 
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-0.04 



< -0.82 < 0.04 

< -0.82 < 0.06 



-+0.18 



-0.70 



+0.24 

0.28 



< 0.06 



TABLE I: Parameter constraints for scalar dark energy models. BASE is WMAP+VSA+CBI+SDSS. All errors are la, upper 
bounds are 2a. 
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< 0.04 

< 0.04 



TABLE II: Parameter constraints for phantom crossing models. BASE is WMAP+VSA+CBI+SDSS. All errors are la, upper 
bounds are 2a. 



7 that controls the importance of the terms involving fi^ 
at late times, 



fid (a) = 



fiS (1 



ni a 3w ° 



fif, 1 



(24) 

When setting 7 = 1, this reduces to our parameteriza- 
tion while increasing 7 shifts the departure from a 
cosmological constant-like behavior of fid (a) to higher 
redshifts. If we require fid (a) to stay monotonic, then 7 
is bounded by 



7 



< 



fid — fijj 
fi|(l-fid)' 



(25) 



Repeating the Monte-Carlo analysis for the scalar dark 
energy case with this extended parameterization (|24D at 
different values of 7 does not alter the constraints on the 
other parameters. It therefore suffices to use the simpler 
parameterization of Equation (jSJ . 
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